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We consider the Fermi quantization of the classical damped harmonic oscillator (dho). In past 
work on the subject, authors double the phase space of the dho in order to close the system at each 
moment in time. For an infinite-dimensional phase space, this method requires one to construct a 
representation of the CAR algebra for each time. We show that unitary dilation of the contraction 
semigroup governing the dynamics of the system is a logical extension of the doubling procedure, 
and it allows one to avoid the mathematical difficulties encountered with the previous method. 



^1" 

o 
o 

(N 

■4— > 

o 
O 



> : 

00 . 

■ 

o ■ 
o : 

■ 

o ■ 

^ : 

» ■ 

G ■ 

s : 
cr 



PACS numbers: 03.65.Yz, 03.65.Fd 

Keywords: damped harmonic oscillator, unitary dilation, open system 



INTRODUCTION 

The damped harmonic oscillator (dho) is a simple clas- 
sical dissipative system which, upon quantization, yields 
a useful example of an open quantum system Un- 
stable particles can be treated as an open system; these 
have been incorporated into quantum field theory in . 
Additionally, modeling a quantum measurement as an 
open system permits one to regard the reduction of a 
state as a continuous process |2|. Both Fermi and Bose 
quantization of the dho have been considered in the lit- 
erature 0, 0, 0, • A technique central to these consid- 
erations is the doubling of the degrees of freedom of the 
classical phase space, allowing one to effectively close the 
system for each moment in time. The extra degrees of 
freedom function as a sink with which the oscillator inter- 
acts and have been related to quantum noise effects 
The doubled system as a whole has been expressed in the 
framework of quantum deformed algebras and related to 
finite-temperature quantum systems 

We shall restrict our discussion to the Fermi quan- 
tized dho, although the techniques are similar in spirit 
for bosons. After introducing the classical dho, we will 
discuss the doubling procedure in terms of the unitary 
dilation of a contraction operator at one instant in time. 
We shall relate this technique to the representation of 
quasifree states over the the canonical anticommutation 
relation (CAR) algebra. Following the model developed 
in , we introduce the unitary dilation of the contraction 
semigroup which describes the time evolution of the dho. 
We then discuss the quantized dho in this framework in 
order to demonstrate the usefulness of this formulation. 



CLASSICAL DHO 

A classical system can be described by a phase space 
M of even (or infinite) dimension whose elements specify 
momenta and position. We denote the complex structure 
on this space by J and the projection onto the momen- 
tum subspace as P. We intend the complex structure 



to generate harmonic oscillations necessitating the addi- 
tional prescription 



JP = (1 - P)J. 



(1) 



For the dho, we take the damping term to be linear in 
the momentum. Damped oscillations arc then generated 
by the operator 



Z = ujJ - 2-yP, 



(2) 



where w is the natural (real) frequency of the oscillator 
and 7 > is the damping strength. The dynamics of a 
particular point in phase space m € M are governed by 
the semigroup 



m(t) = T t m, T t = exp[Zt] 



(3) 



for t > 0. Using the relation in JJJ, we may show that 
the generator satisfies the quadratic 



(Z + 7 1) 2 = ( 7 2 



co 2 )l. 



(4) 



From this identity, we may write the exponential in @ 
in a more tractable form which explicitly demonstrates 
the usual behavior of a dho, depending upon the relative 
values of the damping factor and the natural frequency. 
Defining a = ( 7 2 — uj 2 ) 1 / 2 , we find 



exp[Z< 



cosh ail 



sinh at 



(Z + jl) 



(5) 



Should 7 2 = uj 2 , the case of critical damping, then we 
intend the above relation to be taken in the limit of 
vanishing a; as such, we have 



lim 

a^0 



sinh at 



t. 



(6) 



An underdamped oscillator is characterized by 7 2 < uj 2 
so that a = iu>d is a purely imaginary number. In this 
case, we may write the hyperbolic functions in a more 
transparent form 



coshai = cos to dt, 



sinh at sinu>dt 



(7) 
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hence, one has oscillations of frequency u>d modulated by 
the decaying exponential. For an overdamped oscillator, 
7 2 > uj 2 , one has real a so that the dynamics consist of 
decaying exponentials with two different decay constants. 



CONTRACTION DILATION 

Given Eq. © , it is apparent that T t is a contraction for 
all relevant t and, in fact, strongly converges to zero as t 
tends to infinity. From Ref. @, we extract the following 
result concerning any contraction T on M. First, we 
define the isometric injection j of M into two copies of 
the phase space M = M ®M 



jm = m © 0, 



ni E M. 



(8) 



Then, on the doubled phase space, we may construct 
from the contraction the orthogonal operator 



U 



which satisfies 



T (1 — 7 l 7 l *)i/2 

(X_T*T)V 2 — T* 



j*Uj = T. 



(9) 



(10) 



In words, doubling the dimension of the phase space al- 
lows one to work with an orthogonal operator, instead 
of a contraction. The additional copy of the phase space 
can be regarded as a sink, or coupled oscillator. 

As we are interested in the Fermi jejuantization of the 
dho, we work with the algebra CAR(M) where the phase 
space is considered complex with complex structure J = 
J © — J. The creation operators c(m), linear in their 
argument, satisfy along with their adjoints the CAR 



[c(ifij* ,c(n)] + = (m,n) 



[c(m),c(n)}+=0, (11) 



where [•, •]+ is the anticommutator. By definition, the op- 
erator generates a Bogoliubov transformation of the 
algebra 9]. As this is true for any contraction on M, in 
particular T t , we conclude there is a Bogoliubov transfor- 
mation associated with the dilation Ut relating the Fock 
representation of CAR(M) (at t = 0) and the dho on 
the doubled phase space at any fixed time t. This estab- 
lishes connection with the previous work on the subject 



In general, the operator Ut does not commute with the 
complex structure on the doubled space. We decompose 
the operator into the sum of a complex linear ajj t and 
conjugate linear bu t operators 

a Ut =\{U t - JU t J), b Ut = l(U t + JUj). (12) 

In the usual manner, we define the transformed creation 
operator 



these satisfy the CAR. From [l0|, we find that the Bo- 
goliubov transformation is implementable if and only if 
bu t is Hilbcrt-Schmidt. Given an implementable trans- 
formation, one may construct from elements within the 
Fock representation a vacuum vector that is annihilated 
by Ct(m)* Q; in this sense, the vacuum vector can be 
thought of as a dynamic object. For a finite dimensional 
phase space, this situation is assured. One may show that 
for fixed t > 0, the square of the Hilbert-Schmidt norm 
of bu t scales with the dimension of M. As such, Ut is not 
implementable for an infinite-dimensional phase space. 
As a consequence, one must construct a different repre- 
sentation for each time t with no implementable means 
to change between any two representations. We remedy 
this situation below. 



QUASIFREE STATES 

First, we use similar language to elucidate the con- 
nection between this doubling procedure and general 
quasifree states over CAR(M). A quasifree state ip is 
determined by the two-point correlation functions. The 
state can be characterized by two bounded operators R 
and (conjugate-linear) S. These satisfy < R = R* < 1 
and S* = —S with the two-point correlation functions 
given by 

ip[c(m)* c(n)] = (m,Rn), (p[c(m)c(n)] = (Sm,n), 

(14) 

for m, n £ M . The representation theory of quasifree 
states over the CAR algebra is well developed yjj, 

12]; 

briefly, any quasifree state can be represented as a Fock 
state of the CAR algebra over two copies of M. For the 
most general quasifree state, the connection between its 
representation and the doubling procedure used with the 
dho is most easily accessed through Araki's self-dual rep- 
resentation of the algebra 0. We shall not discuss the 
details of this method because the needed definitions and 
notation would take us too far afield. Instead, we restrict 
the discussion to quasifree states which are invariant un- 
der global U(l) phase changes; these have S = 0, so that 
the state is characterized by R alone. We are justified in 
restricting our scope as these are the relevant states for 
physics. 

Given that, we note from above that R is required to be 
a positive contraction. As such, its square root is defined, 
and \/R is also a contraction. We use the injection j 
defined in set the complex structure J = i © —i, and 
unitarily dilate as before 



c t (m) = c(a Ut m) + c(b Ut m)*; 



(13) 



V = 



Vr VT^r 



(15) 
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With this complex structure, one has 

VR \ , ( 







VT^R 

x/1 - ~R 



(16) 

Using transformed creation operators, as in l|13[) . act- 
ing on the Fock vacuum il, we calculate the two point 
correlation function for elements of M injected into the 
doubled space 

(cv(jm)n, cv(jn)Q) = (c(avjrn)fl, c(avjn)fl) 
= (m, Rn) 

= <p[c{m)*c{n% (17) 

as one would expect. 

The parallel with the quantized dho is particularly ger- 
mane for KMS, or thermal, states. A KMS state ipp, 
at inverse (positive) temperature (3, is a quasifree state 
which describes a quantum system with Hamiltonian H 
in thermal equilibrium [13| . The defining property of this 
state is the commutation relation 



(pp[AB] = (pp[BA 



(18) 



for operators A, B, where we have used a subscript to 
denote time evolution in the Heisenberg picture 



A t = UTAU t 



(19) 



with Ut 



-iHt 



Using the KMS condition {TBI and the 
CAR (|llf) , one can show that the state satisfies the two- 
point correlation functions in l|14f> with R = (l + e - ^) -1 
and S — 0. Hence, we may represent the thermal state 
as a Fock state on two copies of the space for each tem- 
perature (3 > 0, an index reminiscent of time in the dho. 
This connection between the dho and thermal states was 
discussed in 0. 

The Fermi quantized dho and quasifree states share a 
similar mathematical structure with regard to their rep- 
resentation on the doubled space. In particular, the anal- 
ogy between the two is especially compelling for thermal 
states given that they are both indexed by positive num- 
bers (temperature and time). However, the KMS states 
exhibit the much richer structure of Tomita-Takesaki the- 
ory [lij l . As such, they necessarily have a unitary dynam- 
ical component in their definition. When quantizing the 
dho in the above manner, the dynamics of the oscillator 
are, in some sense, frozen out; that is, one has a different 
representation at each moment in time. Heuristically, the 
motivation of the doubling procedure is quite different for 
these two structures. 



SEMIGROUP DILATION 

Above we explicated how one may treat a single con- 
traction as an orthogonal operator. If we are willing to 



inject the original phase space into a space even larger 
than M, then it is possible to unitarily dilate the con- 
traction semigroup T t for all t > 0. In what follows, 
we shall maintain the same notation as above in order 
to make clear that this is a logical extension of the dou- 
bling procedure. The technique that follows can be found 
in |8j ; the application of this theory to the dho was ex- 
pounded upon in The space into which we inject M 
is M = L 2 (R, PM). The injection j : M —> M is given 
by 



(jm)(t) = 2yfyQ(t)PT t m, 



(20) 



with O(t) the Heaviside function. This map can be shown 
to be an isometry. Time translation of the elements of 
the space is given by the unitary operator 



(U t rh)(s) = m{s + t). 



(21) 



Analogous to Eq. I|l()(l . the following holds by construc- 
tion 



fU t j = T t . 



(22) 



Using this technique, one no longer need bother with 
Bogoliubov transformations in the above prescription or 
their implementation. Though M is a much larger space 
than the doubled one, the dynamics of the dho plus en- 
vironment are now unitary. In order to keep track of the 
dho in this space, we employ the projection Q — jj* on 
M. One can then make the orthogonal decomposition 



M = jM (1 - Q)M, 



(23) 



which delineates the oscillator and the environment. This 
method of dilation also exhibits a richer structure than 
that of the doubling procedure. For instance, we may use 
the Fourier transform to ascertain the energy spread of 
an element of jM. As energy is dual to time, we have 



1 



/2tt J- 



-iEt 



(jm)(t)dt. (24) 



After some manipulation [3j , this can be shown to be 



(Fjm)(E) = 



TT E 2 



i2Ej 



P(ujJ + iE)m, (25) 



which has an amplitude reminiscent of the relativistic 
Breit-Wigner amplitude for unstable particles (cf. |l5|). 

One may quantize this total closed system in the usual 
manner by considering CAR(M). The Fock representa- 
tion is now adequate to describe the quantum dho. The 
vacuum vector Q is a stationary state, and the dynamics 
are unitary. We may still address the dissipative nature 
of the dho in this space. For this, we consider the sec- 
ond quantization Q of the projection Q which obeys the 
commutation relation 



[Q, c(m)] = c(Qm) 



(26) 
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for fh £ M. Recalling the time dependence in the Heisen- 
berg picture l|19|) . we calculate the following expectation 
value in the Fock representation 

(c(jm)Cl,Q(t)c(jn)Q) = (jm,Q(t)jn) 

= LfUtjm,j*U t jn) 

= (m,T*T t n), (27) 

for n,m € M. This exhibits the behavior that one would 
expect. 

For the case of the critically damped oscillator, we note 
that given a proper choice of initial conditions one can 
model a collection of unstable particles as in Q • For ele- 
ments of the (nontrivial) subspace N — ker(Z + 7I), the 
evolution of the critically damped oscillator is strictly ex- 
ponential, T t |jv = e -7 *. As a result, for n E N, the time 
dependence of the projection onto the oscillator subspace 
exhibits exponential decay 

(c(jn)n,Q(t)c(jn)n) = e" 2 ^. (28) 

This is another method by which unstable particles can 
be included in quantum field theory. 

CONCLUSION 

In summation, we feel that the method of unitary di- 
lation of the contraction semigroup is a more effective 
means with which to consider the quantized dho. Rather 
than closing the open system at each instant in time by 
doubling the dimension of the phase space, we close the 
entire system at once for all future time. The advan- 
tages of this method are unitary dynamics and the need 
for only one Fock representation with a stationary vac- 
uum vector. The drawback for this approach is that the 
analogy between the dho and thermal states is no longer 
valid; however, we feel that our exposition demonstrates 
that the two systems are conceptually and physically dis- 
parate, lessening the significance of this shortcoming. 
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